Particle-scale reversibility in athermal particulate media below jamming 
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We perform numerical simulations of athermal repulsive frictionless disks and spheres in two 
and three spatial dimensions undergoing cyclic quasi-static simple shear to investigate particle-scale 
reversible motion. We identify three classes of steady-state dynamics as a function of packing 
fraction cj> and maximum strain amplitude per cycle 7max. Point-reversible states, where particles 
do not collide and exactly retrace their intra-cycle trajectories, occur at low <f> and 7 max . Particles 
in loop-reversible states undergo numerous collisions and execute complex trajectories, but return 
to their initial positions at the end of each cycle. Loop-reversible dynamics represents a novel form 
of self-organization that enables reliable preparation of configurations with specified structural and 
mechanical properties over a broad range of <f> from contact percolation to jamming onset at <j)j. 
For sufficiently large <j> and 7max, systems display irreversible dynamics with nonzero self-diffusion. 

PACS numbers: 83.80.Fg61.43.-j,63.50.Lm,62.20.-x 



In equilibrium thermal systems, reversible processes 
are imagined as transitions from one thermodynamic mi- 
crostate to another with no change in free energy. For 
example, thermal fluctuations in high-temperature fluids 
give rise to particle motions that yield only very small 
changes in entropy. However, finite deformations of su- 
percooled liquids and amorphous solids, upon reversal 
of the strain, can produce microscopically irreversible 
motion, such as collective particle rearrangements [![, 
anelasticity, and plastic flow The identification of 
topological defects in crystalline materials is straight- 
forward, whereas it is much more difficult to identify 
particle-scale motion that gives rise to plasticity in amor- 
phous materials 

Granular materials, foams, and other athermal partic- 
ulate media are highly dissipativc, and therefore must be 
driven to induce particle motion. Experimental studies of 
granular media have shown macroscale reversibility of the 
packing fraction during cyclic shear Q and vibration Q . 
Since these systems are far from thermal equilibrium, one 
might assume that they do not display microscale re- 
versible motion when subjected to cyclic driving. Experi- 
mental and computational studies of 2D foams have iden- 
tified both reversible and irreversible TI bubble neigh- 
bor switching events during cyclic shear Q. Researchers 
have also shown that motion of individual particles tran- 
sitions from reversible to irreversible beyond a density- 
dependent critical strain, which decreases with increasing 
packing fraction, in cyclically sheared suspensions 0,0- 

An important open question is whether athermal par- 
ticulate media can undergo completely reversible mo- 
tion due to inter-grain collisions when subjected to cyclic 
loading. We address this question by performing nu- 
merical simulations of frictionless granular materials in 
two and three spatial dimensions undergoing quasistatic 



cyclic simple shear Q over a wide range of packing frac- 
tion <j) and shear strain amplitude 7 max . We identify two 
classes of grain-scale reversible motion, point and loop 
(which are stable to finite perturbations). For point- 
reversible dynamics, particles do not collide during the 
forward cycle, and thus they exactly retrace their trajec- 
tories upon reversal. In contrast, particle collisions occur 
frequently during loop-reversible dynamics, but the sys- 
tem self-organizes so that particles return to the same 
positions at the beginning of each cycle, despite com- 
plex particle motion during the cycle. We map out the 
'dynamical phase diagram' as a function of 4> and 7max- 
The system transitions from point- to loop-reversible and 
then from loop-reversible to irreversible (with nonzero 
self-diffusion) dynamics with increasing <f> and 7 max - We 
show that the time evolution toward steady-state point- 
and loop-reversible behavior can be collapsed onto a uni- 
versal scaling function with power-law scaling at short 
and intermediate times, and exponential decay at long 
times. 

We perform numerical studies of N athermal disks and 
spheres undergoing quasi-static, cyclic simple shear in 
2D and 3D at constant <p using shear-periodic boundary 
conditions in square or cubic cells (Tlj . Particles interact 
via the pairwise, purely repulsive linear spring potential 

n^) = |(i-^|) e(^-^), (i) 

where is the center-to-center separation between par- 
ticles i and j, Q(x) is the Heaviside step function, <7y = 
((7; + <Jj)/2, and <Ji is the diameter of particle i. We fo- 
cus on bidisperse particle size distributions, i.e. 50-50 
mixtures by number with diameter ratio <ji/<j s = 1-4, to 
frustrate crystallization during shear [l2|- We consider 
system sizes from TV = 32 to 512 to assess finite size ef- 
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fects for packing fractions below and near the onset of 
jamming (</>,/ - 0.84 [n| in 2D and - 0.65 [U] in 3D). 

The particles are initially placed randomly in the sim- 
ulation cell at packing fraction <f> and then relaxed using 
conjugate gradient energy minimization (l2T |. We apply 
simple shear strain by shifting each particle horizontally 
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in increments of A7 = 10~ 3 , where 
coordinates of particle i at the fcth step of strain cycle 
n After each strain step, we minimize the total 

potential energy at fixed shear strain, jk = kAj for the 
forward or 7^ = 27, nax — kAj for the reverse part of the 
cycle. This process is repeated for up to n = 10 6 cycles. 

During the simulations, we measure the single cycle 
mean-square displacement (at step k = for each cycle) 
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are the non-affine displacements of parti- 
cle i after subtracting off the affine contribution. The 
long-time dynamics arc either reversible or irreversible 
depending on <f> and 7 max . We quantify the steady-state 
behavior by measuring Ar 2 (n) and L 2 (n) of the intra- 
cycle particle trajectories. 

Particles in point-reversible systems organize to avoid 
collisions. At long times, no collisions take place, and 
L(n) = Ari(n) = 0, or more aptly, they fall below small 
numerical thresholds, e.g. Ari(n) < t t = 5 x 10~ 4 and 
L(n) < tl = 10~ 8 . The values of t t and tl do not qual- 
itatively affect our results as long as they are sufficiently 
small. Particle motions for point-reversible systems are 
affine and in the direction of the imposed affine shear 
(Fig. [JJ (a)). Thus, the non-affine tracks of each particle 
are zero (Fig. [2(d)). 

In loop-reversible systems, particle collisions occur fre- 
quently, but the system self-organizes so that particles 
return to the same positions at the start of each cy- 
cle. Since collisions between particles occur, L(n) > 0, 
but Ari(n) = (i.e. below r r ). Particle trajectories 
(Fig. Q] (b)) and non-affine displacements (Fig. [T] (e)) 
trace out complex paths, yet all particles end up in the 
same locations at the beginning of each new cycle, i.e. 
X 



n+1,0 



Y 7 l in 2D. Thus, particle 
trajectories form closed loops in configuration space. We 
focus on period one loop-reversible systems, but multi- 
period dynamics are also found. Particle trajectories are 




FIG. 1: (color online) Particle tracks (solid lines) in a small 
window of a N — 128 system of bidisperse disks undergoing 
(a) point-reversible, (b) loop-reversible, and (c) irreversible 
behavior during cyclic simple shear. Panels (d), (e), and (f) 
show the disks' tracks in panels (a), (b), and (c), respectively, 
after subtracting off the affine motion. The systems in (a) 
and (d) correspond to (j> = 0.64 and 7 max = 0.8, (b) and (e) 
to (f> = 0.8 and 7 max = 0.8, (c) and (f) to (j> = 0.82 and 
7max = 0.8. Pluses, circles, and crosses mark the beginning, 
middle, and end of the particle tracks, and particle outlines 
correspond to the beginning of the cycle. 



elongated in the direction of affine displacement, whereas 
nonaffinc displacements are more compact. 

Particles in systems undergoing irreversible dynamics 
do not return to their original positions at the beginning 
of each new cycle, e.g. X l n+10 ^ X l n and Y* +10 ^ 
Y^q in 2D. (Figs. [T] (c) and (f).) Irreversible systems 
have nonzero Ari(n) and L(n) (i.e. L(n) > T£ and 
Ari(n) > r r ). Systems can be "transient" irreversible 
in time and evolve into point- or loop-reversible systems, 
or steady-state irreversible and remain irreversible in the 
large-cycle limit with nonzero self-diffusion. 

The steady-state "dynamical phase diagram" in Fig. [2] 
(a) for cyclically sheared athermal disks shows point- 
and loop-reversible , as well as irreversible regimes ver- 
sus <fi and 7 m ax [14| . Point-reversible systems occur at 
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FIG. 2: (color online) (a) "Dynamical phase diagram" 
for N = 128 bidisperse disks after n = 10 4 cycles show- 
ing point-reversible (circles), period-one (pluses) or multi- 
period (crosses) loop-reversible, and transient (triangles) or 
steady-state (squares) irreversible dynamics versus j mBX and 
<j>. The solid and dashed lines indicate <^>L(7max) and </>i(7 m ax), 
the boundaries between point- and loop-reversible dynamics 
and between loop-reversible and irreversible dynamics, re- 
spectively. The vertical dotted and dot-dashed lines define 
4>r — 0.785 and (j>j « 0.84. (b) Intra-cycle mean-square 
displacement Ari versus arc-length L after n = 10 (blue), 
3 x 10 2 (red), and 10 4 (black) cycles. The regions labeled 
P, L, T, and / define point-reversible, loop-reversible, tran- 
sient irreversible, and steady-state irreversible dynamics, re- 



spectively. The P-region extends from L = A 



n 



10" 



to 



10 _ (not shown). The dashed boundaries indicate An = r r , 
An = 0.3cr s , L — tl, L — 30Ari, and L — An discussed 
in the main text. The inset shows the fraction Ft of systems 
in (a) categorized as point-reversible (solid), loop-reversible 
(dashed), transient irreversible (dotted), and irreversible (dot- 
dashed) dynamics versus n. Similar data for 3D systems is 
shown in Supplementary Material. 



low <p an d 7max, whereas irreversible systems occur for 
</> > <j)j. At intermediate packing fractions between 
roughly contact percolation [10[ and jamming onset, e.g. 
0-6 % <t> ^$ 0.84 in 2D, loop-reversible systems are found. 
The boundary between point- and loop-reversible sys- 
tems is 7 max ~ A((j))((j)j - cj)) x e{(j)j - 0), where A(4>) 
depends weakly on cf> and A ~ 1.2 ± 0.1 for <f> — > cf>j 



and 2.2 ± 0.2 for c/> <C <fij. Over a finite number of cy- 
cles (i.e. n < 10 4 in Fig. [2] (b)), transient irreversible 
dynamics can occur, but these systems become point- 
reversible, loop-reversible, or steady-state irreversible as 
n — > oo. Point-reversible systems tend to form ordered, 
size-segregated layers, in which particles cannot collide 
during simple shear. Further, the loop-reversible to irre- 
versible transition in steady-state <^_R(7max) is bounded 
in the large-7 max limit by the highest packing fraction 
at which systems can form ordered, size-segregated 
layers; </>,r(7 max — > oo) = 4>r = 7r/4 ~ 0.785 in 2D and 
0.605 in 3D in the N -> oo limit. 

In Fig.[2](b), scatter plots of L(n) versus Ari(n) for 2D 
systems illustrate the evolution of the dynamics with in- 
creasing n. The points form several well-defined clusters: 
point-reversible (P) with L < tl and Ari < r r , loop- 
reversible (L) with nonzero L (L > tl) and Ari < r r , 
and irreversible (/) with nonzero L (L > tl) and Ari 
(Ari > T r ). The P, L, and I clusters are separated by 
more than 3 orders of magnitude in Ari or L. For re- 
gion L, we also mandate L > 30Ari since systems with 
L < 30Ari typically relax to point-reversible states. We 
also enforce Ari > 0.3a s to define region I since sys- 
tems with Ari < 0.3c s typically relax to point- or loop- 
reversible states. Systems that do not fall within the 
bounds defining regions P, L, and I are categorized as 
transient irreversible (T). As n increases, the fraction F t 
of systems in the transient regime vanishes as a power-law 
n~ a (where a — 0.56 ± 0.01), while the fraction of point- 
reversible, loop-reversible, and steady-state irreversible 
systems saturates near 10 4 cycles (inset to Fig.[2](b)). 

We also characterized the dynamics of these systems 
as they approach steady-state point- and loop-reversible 
states (Fig. [3] (a)). We find that the single cycle mean- 
square displacement can be described by a function that 
interpolates between power-law and exponential decays 
at short and long times, respectively: 



Ar{>) = f+(n)(n/n c )- a + f_(n)e 



(5) 



where f±(n) = (1 + e ±7(n-n c )^-i ^ ^ ^ ^ Bc j s ^ c C y C i c 
number at which the decay changes from power-law to 
exponential behavior, a is a power-law scaling exponent, 
and (3 characterizes the exponential decay. 

In Fig. [3] (b), we plot the best fit Ar{(n) versus 
Ari(n) at each 7 max and <fi (averaged over 16 initial 
conditions) for the 2D systems in Fig. [2] that evolve to 
point-reversible states. The scaling function in Eq. [5] col- 
lapses more than 60% of point-reversibile systems with 
deviations A = ((log 10 Ar( (n) - log 10 An(n)) 2 ) < 0.18. 
The top and bottom insets in Fig. [3] (b) show the 
power-law scaling and exponential decay of Ari (n) sep- 
arately. We find similar scaling for the approach to loop- 
reversible states. However, the exponential decay for 
loop-reversible systems is difficult to differentiate from 
numerical error because the long-time dynamics occurs at 
larger n c and smaller Ari than that for point-reversible 
systems. In Fig. [3] (c), we show the power-law decay 
for all systems that evolve to loop-reversible dynamical 
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a log. 



FIG. 3: (color online) (a) Single cycle mean-square displace- 
ment Ari versus n for a point-reversible system at <f> = 0.64 
and 7 m ax = 0.5 (circles) and a loop-reversible system at 
4> — 0.8 and 7 max = 0.5 (pluses) with best fits to Ar[ 
(Eq. [5]) indicated by solid and dashed lines, (b) Compari- 
son of An(n) (averaged over 16 initial conditions for each 
7max and (f>) to Ar{(n) (black dots), where n is the cycle 
number, for point-reversible systems in Fig.[2]with A < 0.18. 
Ar((n) = Ari(n) is indicated by the dashed line. The top left 
inset shows log 10 Ari (n)/o versus a log 10 n/n c (black dots). 
The dashed line indicates Ari(n)/a = (n/n c )~ a . The bot- 
tom right inset shows log 10 Ar\{n)/a versus fi{n — n c ) (black 



dots). Ar\(n)/a 



-/3(n-n c ) 



is indicated by the dashed line. 



(c) log 10 Ari in) / a versus a log 10 n (black dots) for systems 
in Fig. [5] that evolve to loop-reversible dynamical states with 
A < 0.04. Ari(n)/a — n~ a is indicated by the dashed line. 
The inset shows An /a versus n for three independent initial 
conditions (squares, triangles, and pluses) at (j> = 0.76 and 
7max = 0.8. Exponential fits to the large n regime are shown 
as solid, dashed, and dotted lines with slopes f3 = 0.029, 0.026, 
and 0.016, respectively. 



states. In the inset, we also show several systems for 
which wc captured the long-time exponential decay. 

In Fig. [4] (a), we show the power- law scaling exponent 
a for 2D systems that evolve to point- and loop-reversible 
states versus <j> and 7 max - Wc find that a < 1 for all loop- 
reversible systems and point-reversible systems near the 
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FIG. 4: (color online) (a) Contour plot of the power-law expo- 
nent q (Eq.[S]) versus <j> and 7 max for 2D systems that evolve to 
either point- (circles) or loop-reversible (crosses) states, (b) 
Contour plot of the cycle number n c (Eq. [5]) that controls the 
crossover from power-law to exponential decay versus cj> and 
7max for the point-reversible systems in (a). 



crossover from point- to loop-reversible behavior, which 
suggests that the origin of the slow dynamics is related to 
contact or "collision" percolation. In contrast, a > 1 for 
point-reversible systems at low <j> and 7 max . In Fig.@](b), 
we plot n c for point-reversible systems versus <\> and 7 max . 
We find that n c increases with <f> and 7maxi an d appears 
to be diverging as the system approaches the transition 
from point- to loop-reversibility. 

To test the stability of steady-state loop-reversible 
states, we perturb all particles at strain 7 = by an 
amplitude 5 in random directions [l4| . We then perform 
cyclic simple shear on the perturbed system and measure 

the deviation, A r = J (Na^ 



where 



Si Vnfl ' n,0| 

' n are the coordinates of the perturbed system after 
t cycles required to reach steady steady at each <j> an d 
7max- We find that loop- reversible systems are stable 
(A r < T r ) for perturbations 8 < S c ~ 10" , i.e. per- 
turbations as large as one-tenth of a particle diameter, 
where 6 C is relatively insensitive to cf> for 7 max ^ F 

In conclusion, we studied the extent to which particle- 
scale motion can be reversible in athermal systems under- 
going cyclic quasistatic loading. Wc identified two types 
of reversible behavior. For point-reversible states, par- 
ticles do not collide and therefore trivially retrace their 
paths. For loop-reversible states, all particles undergo 
multiple collisions and have complex trajectories, yet all 
particles return to the locations they were in at the begin- 
ning of the cycle. Wc determined the regions of packing 
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fraction and strain amplitudes in 2D and 3D where these 
dynamical states are stable. In particular, we find that 
loop-reversible states occur over a range of packing frac- 
tion from contact percolation [Io| to jamming onset, and 
thus our results emphasize that complex spatiotemporal 
dynamics are found well below <f>j. Loop-reversible dy- 
namical states represent a novel form of self-organization 
that will enable reliable preparation of configurations 



with particular structural and rheological properties over 
a broad range of packing fractions. 
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